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We present a deduction of the Lindblad equation that is accessible to students with a basic 
background on quantum mechanics. We consider a specific case, corresponding to a very simple 
situation, where a primary system interacts with a bath of harmonic oscillators at zero temperature, 
with an interaction Hamiltonian that resembles the Jaynes-Cummings format. We start with the 
Born-Markov equation and, tracing out the bath degrees of freedom, we obtain an equation in the 
Lindblad form. The specific situation is very instructive, for it makes it easy to realize that the 
Lindblads represent the effect on the main system caused by the interaction with the bath, and 
that the Markov approximation is a fundamental condition for the emergence of the Lindbladian 
operator. The formal derivation of the Lindblad equation for a more general case requires broader 
considerations regarding the environment, the temperature, and the use of quantum dynamical 
semi-groups than we have considered in the particular case treated here. 



I. INTRODUCTION 

The Lindblad equation [Q is the most general form for a Markovian master equation, and it is very important for the 
treatment of irreversible processes, from dissipation and decoherence to the quantum measurement process |3, 
Contrasting against its importance and wide range of applications, its original deduction |^ involves the forma ism 
of quantum dynamical semigroups ^, which is quite unfamiliar to most of the students and researchers. Other more 
recent ways to derive it involve the use of Ito stochastic calculus |7[ or, in the specific case of quantum measurements, 
considerations about the interaction between the system and the meter These methods, their assumptions, their 
applications, and, more importantly, their physical meanings appear very intimidating to beginning students. 

To make the Lindblad equation more understandable, this article presents its deduction in the specific case of 
two systems: S , the principal system, and B, which can be the environment or the measurement apparatus, at zero 
temperature, with an interaction between them that resembles the one of the Jaynes-Cummings modelH]. Initially 
we derive the Born-Markov master equation |2] and then we trace out the degrees of freedom of system B. The 
Lindbladian emerges naturally as a consequence of the Markov approximation. Each Lindblad represents the effect 
on system S caused by the S — B interaction. 

Clearly, the present approach does not prove the general validity of the Lindblad equation. Our intention is simply 
to provide an accessible illustration of the validity of the Lindblad equation to non-specialists. The only prerequisite 
to follow the arguments exposed here is a basic knowledge of quantum mechanics, including a familiarity with the 
concepts of the density operator and the Liouville-von Neumann equation, at the level of Ref. |l0[| , for example. 

The paper is structured as follows: in the Sec. II we derive the Born-Markov master equation by tracing out the 
degrees of freedom of system B, starting from the Liouville-von Neumann equation; in Sec. Ill we derive the Lindblad 
equation; and in Sec. IV we present the conclusion. 



II. THE BORN-MARKOV MASTER EQUATION 



Let us consider a physical situation where a principal system S, whose dynamics is the object of interest, is coupled 
with another quantum system B^ called bath . Here, T-Ls and T-Lb are, respectively, the Hilbert spaces of principal 
system S_ and bath^; the global Hilbert space S + B will be represented by the tensor-product space Hs ® Hb- The 
total Hamiltonian is 



H{t)=Hs(^iB + is®HB+aHsB, (1) 

where Hs describes the principal system S, Hb describes the bath B, Hsb is the Hamiltonian for the system-bath 
interaction and 1b and I5 are the corresponding identities in the Hilbert spaces. For the sake of simplicity, let us 
ignore the symbol ® and write 



H{t) = Hs + HB + aHsB. 



(2) 
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Here, a is a real constant that provides the intensity of interaction between the principal system and the bath. Writing 
PSB for the global density operator {S + B), the Liouville-von Neumann equation will be: 



Hs + Hb + aHsB , PsB 



(3) 



It is convenient to write Eq. in the interaction picture of Hs + Hb- With the definitions of the new density 
operator and Hamiltonian: 



and 



the new equation forp (t) will be 



H{t) = e^i"'+"^yHsBe-^^"'+"^> 



p{t) = eU"'+"-ypsB it) e-U"s+HB)t^ 



^^p{t) = ~'-a Hit), Pit) 



We want to find the evolution for pg it) ^ trB {psB where 



(4) 



(5) 



(6) 



PSB it) = e-ii^^+"-ypei(^-+"->. 
Equation (^) is the starting point of our iterative approach. Its time derivative yields 



(7) 



pW = /5(0)-^« I [Hit'), Pit') 



dt'. 



Replacing Eq. (^) into Eq. (||), we have 



dtP^'^ = -r 



Hit),piO) 



ft2 



Hit), Hit'), Pit') 



dt' 



(8) 



(9) 



For th e Born approximation , Eq. (^) is enough. Then, we take the partial trace of the bath degrees of freedom, 



^/5s it) = -^'^trB { [h it) , p (0)] } - ^ahrB | ^ (0 , ^ 



Hit), I Hit'), Pit') 
'- 



dt' 



(10) 



By the definition in Eq. H (t) depends on Hsb, and Hsb can always be defined in a manner in which the first 
term of the right hand side of Eq. (foh is zero. Hence, we obtain 



d . , , 1 2 
JtPsit)--j^^ trB 



t ^ 



Hit), I Hit'), Pit') 





dt' 



(11) 



Integrating Eq. (O) from t to t' yields 



PS it) -PS it') = --^a' / dt'tvB 



Hit'), Hit"), Pit") 



dt" 



which shows that the difference between ps it) and ps it') is of the second order of magnitude in a and, therefore, 
we can write ps (t) in the integrand of Eq. ([ll[), obtaining a time- local equation for the density operator, without 
violating the Born approximation: 



d ^ , , 1 2 

= trB 



Hit), Hit'), Pit) 
Jo 



dt' 



(12) 
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The a constant was introduced in Eq. (||) only for clarifying the order of magnitude of each term in the iteration 
and, now, it can be supressed, that is, let us take a = 1. Thus, let us write: 



H(t), H{t'),p{t) 
Jo '- 



dt' 



(13) 



For this approximatin, we can write p (t) = ps (t) (E) Pb inside the integral and obtain the equation that will be used 
in the next calculations (again, for the sake of simplicity, let us ignore the symbole ®): 



dt 



dt' 



'trB[[H {i),[H [t'),ps{t)pB]]] 



(14) 



where we assume that the integration can be extended to infinity without changing its result. Equation (|lj) is the 
Born-Markov master equation H]. 



III. LINDBLAD EQUATION 
A. The master equation commutator 

Now, we consider that S commutes with i/g, i. e.. 



S.Hs 



0, 



resulting in 



S{t) ^ S 



(15) 



[S is not affected by the interaction-picture transformation). We will keep S and B [t) as a general operator. The 
later can be defined as: 



B{t) = eJ^^^Be-'*" 
With this, the commutator in Eq. (|l|), H {t) , H (f ) , ps (t) pB 

H{t),\H {t'),pBPs (t) 



will be evaluated. Firstly, 



= h 

= h 
+ h 



SB^ {t) + S^B{t),[H {t'),pBPs (t) 
SB^ it),[H {t'),pBPs (t) 
S^B{t),\H {t'),pBPs (t) 



The gradual expansion of each term in Eq. (|17|) will result in 



(16) 



(17) 



SB^t), H{t'),pBPs{t) 



= h 



SB^ [t] , \SB^ (f) + S^B it') , pBPs (t) 



hSB^ it) SB^ [t') + S^B it') , PBPS (t) 



SBHt') + S^B{t'),pBPs{t) SB^t) 



hSB^ (t) \SB^ (f) + S^B {t')] PBPS (t) 
hSB^ (t) PBPS it) \SB^ (f) + S^B (t') 



SBHt') + S^B{t') pBPs{t)SBHt) 



hpBPs (t) \SB^ it') + S^B (t')] SB^ (t) 



(18) 



and 
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S^B (t) , [h it') , pBPs (t)] ] = h [S^B (t) , [SB^ (f) + S^B [t') , pbRs (t) 

= hS^B (t) [SB^ (f) + S^B it') , PBPS (t) 

- h \SB^ it') + S^B (t') , PBPS (t)] S^B (t) 
= hS^B (t) [SB^ (t') + S^B (t')] PBPS (t) 
^ hS^B (t) PBPS (t) [SB^ (<') + S^B (t') 

- h \SB^ it') + S^B (t')] PBPS (t) S^B (t) 
+ hpBPs [t) [SB^ it') + S^B (t')] S^B it) , 

or, expanding Eqs. (|l8|) and ( pj| ) and grouping the similar terms in S and B, we have: 



(19) 



SBHt), Hit'), PBPS (t) 



hSB^ it) SB^ it') PBPS it) - 
hSB^ it) PBPS it) SB^ it') - 
hSB^ it') PBPS (t) SB^ it) 
hpBPs it) SB^ it') SB^ it) - 
hSSps it) St it) St it') pB 
hSps it) SB^ it) pbB^ it') - 
hSps it) SB^ it') pbB^ it) 
hps it) SSpbB^ it')B^ it)- 



^ hSB^ it) S^B it') PBPS it) 

- hSB^ it) PBPS it) S^B it') 
hS^ Bit') PBPS it) SB^ it) 

- hpBPs it) S^B it') SB^ it) 
~ hSS^Ps it) St it) B it') PB 

- hSps it) S^B^ it) pbB it') 

- hS^ps it) SB it') pbB^ it) 
-hps it) S^SpBBit')B^ it) 



(20) 



and 



S^Bit), Hit'), PBPS it) 



hS^Bit) SB^ it') PBPS (i)4 
hS^Bit) PBPS it) SB^ it')- 
hSB^ it') PBPS it) S^Bit)- 
hpBPs it) SB^ it')S^Bit)^ 
hS^Sps it) Bit) B^ it')pB- 
hS^ps it) SBit)pBB^ it')- 
h^Sps it) S^B^ it')pBBit) 
hps it) SS^pbB^ it') B it) 4 



hS^Bit) Bit') PBPS it) 

hS^B it) PBPS it) S^B it') 
hS^B it') PBPS it) S^B it) 
hpBPs it) S^B it') S^B it) 

- hS^S^ps it) B it) B it') PB 
hS^Ps it) S^B it) pbB it') 

- hS^ps it) S^B it') pbB it) 
hps it) S^S^pBBit')Bit). 



(21) 



B. The partial trace 

Now we are in a position to trace out the bath degrees of freedom in Eqs. (^0|) and (^l[). Here, since we do not 
specify the bath and its Hamiltonian, we make a very simple and reasonable hypothesis: we suppose that 

trB[Bit)Bit')pB] = trB[B^ it)B^ it')pB] =Q,^t,t'. 



As wc will show bellow, this is the case if we suppose a bath of bosons. 
With this, then. 



and 



tvB { [SB^ {t) , [h {t') , pBPs (t)] ] } = hSS^ps (t) tre 

- hSps{t)ShrB 

- hS^ps{t)StrB 
+ hps{t)S^StrB 

tTB { [S^B (t) , [h {t') , PBPS {t)] ] } = hS^Sps (t) trB 

- hS^ps{t)StrB 

- hSps{t)ShrB 
+ hps{t)SShrB 



fit {t)B{t')pB 

{t) pbB (f) 
B{t')pBB^ {t) 
PbB (f) St (t) 

B (t) St (t') pB 
B{t) pbB^ it') 

{t')pBB{t) 

PbB^ {t') B (t) 



where, if we use the cichc properties of the trace, 



trB{[sB^ {t),[H{t'),pBPs (i)]]} 



and 



tTB {[s^B{t),[H{t'), PBPS (t)]]} 



h 
h 

h 
h 



SS^Ps (t) ~ S^ps (t) S] tTB {st (t) B (f) pb} 
PS (t) S^S - Sps (t) 5t] tTB {b (f) St (t) pb} 

S^Sps (t) - Sps (t) 5t] trB {b (t) {t') pb] 
PS (t) SS^ - S^ps (t) S] trB {b^ {t') B {t) pb} . 



The terms represented by Eqs. ( p^ ) and ( |25| ) aUow us to return to the Eq. (p^): 



trB{[H{t),[H{t'),pBPs it)]]} 



trB{[sB^ {t),[H{t'),pBPs (t)]]} 
Bit), [nit'), PBPS it)]]} 
SS^Ps it) - S^ps it) S] trB {b^ it) B it') pb} 
PS it) S^S - Sps it) S^] trB {b it') B^ it) pb} 
S^Sps it) - Sps it) S^] trB {b it) St {t') pb} 
PS it) SS^ - S^ps it) S] trB {b^ it') B it) pb} 



hhrB { 



C. The expansion of the integrand of the master equation 

With the rcsuhs of the preceding paragraphs, the integrand in Eq. ( p^ becomes: 

trB [[nit), [nit'), PBPS it)]]} = 



+ 
+ 
+ 



SS^Ps it) - S^ps it) S] trB [b^ it) B it') pb} 
PS it) S^S - Sps it) 5t] trB [b it') B^ it) pb} 
S^Sps it) - Sps it) S^] trB {b it) St {t') pb} 
PS it) SS^ - S^ps it) S] trB {b^ it') B it) pb} 
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For convenience, let us define the functions 

F{t) 
G{t) 



Then, 



dt'trB {b (t) {t')pB}, 
'trB{B^ {t')B{t)pB}. 



F* (t) 
G* it) 



df 



dt 



(28) 



dt 



'trB{B{t')B^ {t)pB} 
'trB{B^ {t)B{t')pB} 



Replacing Eq. (^) in Eq. (|T|) yields: 



SS^Ps (t) - S^ps (t) S\ a* (t) ~ [ps (t) S^S ~ Sps (t) F* (t) 
S^Sps (t) - Sps (t) S^] F (t) - \ps (t) SS^ - S^ps (t) S] G (t) . 



(29) 



Actually, the usual Lindblad equation emerges when G{t) = and F{t) = F*{t). In the following, we make some 
specifications about the environment to discuss these approximations in detail. 



Let us consider the bath hamiltonian 



and the system-bath interaction 



D. The bath specification 



Hb = ti^iOka\dk 

k 



HsB = h (SB^ + S^b) , 



(30) 



(31) 

where dk e are the annihilation and creation bath operators, the Uk are the characteristic frequencies of each mode, 
S" is a general operator that acts only on the principal system S', and B is an operator that acts only on the bath B, 



B 



^Qk^k, 

k 



where gk are complex coefficients representing coupling constantes. Let us consider the interaction picture, 

B{t) = e^^«*i?e~^^^*. 
Expanding each exponential and using the commutator relations, Eq. (|3^ ) will result in 



B{t)^Y.9ldke 



Furthermore, for the initial state of the thermal bath, we consider the vacuum state: 

PB = (|0)|0)...)®((0|(0|...). 



(32) 



(33) 



(34) 



(35) 



The evaluation of the F (t) and G (t) functions defined in Eq. ( |2^ ) are done considering the B (t) and pB definitions 
in Eqs. (|4|) and (||). By Eq. (g), i?t {t) is 

^t(t) = ^5,ate-'=*. (36) 

k 

Then, the partial trace in F (t) and G (t) can be evaluated: 

trB {b it) St (t') pb} = trB {b (t) B^ {t') (|0) |0) ...) ® ((0| (0| ...)} (37) 

and 

trB {b^ it') B (t) pb} = trB {b^ (f) B (t) (|0) |0) ...) ® ((0| (0| ...)} • (38) 
If we use some bath state basis{|6)}, Eqs. ( ^ and ( |38|) become 

trB{B{t)BHt')pB} = ^(6|S(t)Bt(t')(|0)|0)...)®((0|(0|...)|6) 

b 

= Y.m{o\...)\b){b\B{t)BHt')m\o)...) 

b 

= m{o\...)Y,\b){b\B{t)BHt')m\o)...) 

b 

= m{0\...)B{t)BHt')m\0)...) (39) 

and 

trB {sHt') B (t) pb} = ^(&|snO^(t)(|0)|0)...)®((0|(0|...)|6) 

b 

= ^((0|(0|...)|6)(6|Bt(t')B(t)(|0)|0)...) 

b 

= ((0|(0|...)^|6)(6|Bt(t')B(t)(|0)|0)...) 

b 

= m{0\...)BHt')Bit){\0)\0)...). (40) 
Let us expand B'' (t) and B (t) using Eqs. (|4|) and (||): 

irs {S it) fit (f) p^} = ((0| (0| ...) ^ 5^afce— ^ g,.aI,e-'='*' (|0) |0) ...) 

k k' 

= E ff^^fe'e-^"'-*-'^'"*') ((0| (0| ...) a,al, (|0) |0) ...) (41) 

k,k' 

and 

trB {st (t') B (t) pb} = ((0| (0| ...) E 5..al.,e-^'*' ^ .g^a^e— ^* (|0) |0) ...) 

k' k 

= E ((0| (0| ...) 4,a, (|0) |0) ...) = (42) 

k,k' 

Hence, we can rewrite Eq. ( |4l| ) with the a| operators on the left of the operators. We know that 



aka\, = 4,fc' + a\,ak- (43) 

Then 

tvB {b {t) St {t') ps] = Y.9l9k'e-<--'-^''')5k,k' + ^ glgt'e^-"'--''') ((0| (0| ...) at,a, (|0) |0) ...) 

k,k' k,k' 

= 5]|gfc|'e-'"'-(*-*'). (44) 

k 

Therefore, from Eqs. (|2|) and (|^, it follows that 

Fit) = Eiffel' t dt'e-^'-'-i'-''), 
k -^0 

G'(i) = 0. (45) 

E. Transition to the continuum 

In the expression of F (t) in Eq. ( p5[ ) , if we adopt the general definition of the density of states as 

J{uj) = Y,\9i\^H^-^i), (46) 
then the sum over k can be replaced by an integral over a continuum of frequencies: 



F 



{t) = r dujJ{iu) f dt'e-"^ 
Jo Jo 



Let us introduce the new variable 



with 



T = t-t', 
dr — —dt'. 



t pO pt 

dt' = ~ dr = dr. 



yielding 



Jt JO 



/•OO /'t 

F(t) ^ / dujJ{Lu) dTC-"^''. 
Jo Jo 



F. The Markov approximation 



In the Markov approximation, the limit i — > oo is taken in the time integral, as we have mentioned regarding Eq. 
(|l^), that is, we take 
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As the integrand oscilates, we will use the device: 



dre-"^^ = lim / dre^'"^""^ 
'(-^0+ Jo 

1 



lim 
= lim 
lim 



r)^0+ T] + iuj 



i)-i-o+ rf + 



— lim 



7r5 (oj) — iP- 



where P stands for the Cauchy principal part. Then, 



F = TT dujj (w) 5 {u) -iF duj 



G. The final form 



For a general density of states, F yields 



F 



J + is 



where 



7 = 27r / dujJ (cj) S {lu) , 
Jo 

s ^ -2P rdc.^. 

Jo 



As we have verified that G = 0, let us replace Eq. (E^) in Eq. (C 



d 



PS it) = - ps{t)S^S~Sps{t)S^ 



7 — «e 



S^Sps {t)~Sps it) 



7 + «e 



2 



PS (t) S^S - Sps (t) + S^Sps (t) - Sps (t) 
PS (t) S^S - Sps (t) - S^Sps (t) - Sps (t) 5t 



(47) 

(48) 
(49) 



If the density of states is chosen to yield £ = (a Lorentzian , for example, where we can extend the lower limit of 
integration to —00), the final result is: 



dt 



PS (t) = 7 



Sps {t)S^ ~^[s^S,ps (t)} 



Let us, then, return to the original picture. Since 

PS (t) 

then 



(50) 



(51) 



dt 



ps{t) = '-ei^stH^p^^-ifisi + ^iHst^^-i 



Hst 



Hst'iPS 

dt 



e R 



Hst 



Hs,ps 



e ft 



h 

Hst 



eT^"^*psHse~^"''* 



(52) 
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Performing the same operation on the right-hand side of Eq. (pw gives 



Replacing Eqs. 



Sps (t) S^^^[s^S,ps (t)} 
and (pSh in Eq. (pQ), we obtain: 



Hst 



SpsS^ - ^ 15^5,^5! 



e R 



Hst 



(53) 



dps _ _ i_ 
dt ~ h 



Hs,ps 



+ 7 



SpsS'' - i js'^'^S', psj 



(54) 



IV. CONCLUSION 



In summary, in this paper we consider an interaction that resembles the Jaynes-Cummings interaction [Q] , Eq. ( |31 
between a bath and a system S, assuming that the operator S commutes with the system Hamiltonian, Hs, at zero 
temperature. We substituted Eq. (^ij) in the Born-Markov master equation (|lj) and took the partial trace of the 
degrees of freedom of B. The T = hypothesis is necessary to simplify the calculations, making them more accessible 
to the students, simplifying the treatment of Eqs. (B^ and (pq), and avoiding complications such as the Lamb shift 



in Eq. (47). The Markov approximation in Sec. III-F was also vital to obtain the final result, Eq. (p4|). All these 
simplifications limit the validity of our derivation to more general cases, but it provides a detailed illustration of the 
physical meaning of each term appearing in the Lindblad equation. 

Equation (|5j) is commonly presented with several S operators, usually denoted by L, in a linear combination of 
Lindbladian operators. The L operators are named Lindblad operators and, in the general case, the Lindblad equation 
takes the form: 



dps 
dt 



H: 



s,Ps 



Y.\l^PsL]-1{l'A^Ps} 



(55) 



The Lindblad operators represent the effects on system S caused by the S — B interaction, Eq. ( pi] ) . It is important 
to emphasize that, due to our simplifying assumptions, the summation appearing in Eq. ( p5|) was not obtained in our 
derivation of Eq. (^4|). If we consider only the first term on the right hand side of Eq. (p5|), we obtain the Liouville-von 
Neumann equation. This term is the Liouvillian and describes the unitary evolution of the density operator. The 
second term on the right hand side of the Eq. (|54|) is the Lindbladian and it emerges when we take the partial trace - 
a non-unitary operation - of the degrees of freedom of system B. The Lindbladian describes the non-unitary evolution 
of the density operator. 
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